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Summary : We revisit the problem of mutually unbiased measurements in the context 
of estimating the unknown state of a rf-level quantum system, first studied by W. K. 
Wootters and B. D. fields[7] in 1989 and later investigated by S. Bandyopadhyay et al 
[3] in 2001 and A. O. Pittenger and M. H. Rubin [6] in 2003. Our approach is based 
directly on the Weyl operators in the L 2 -space over a finite field when d = p r is the 
power of a prime. When d is not a prime power we sacrifice a bit of optimality and 
construct a recovery operator for reconstructing the unknown state from the probabilities 
of elementary events in different measurements. 
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1 Introduction 

This is almost an expository account of a well-known problem of quantum probability 
and statistics arising in the context of quantum information theory. There is a c?-level 
quantum system whose pure states are described by unit vectors in a <i-dimensional 
complex Hilbert space Ti equipped with the scalar product (<p\ip) between elements tp, i/j 



in H. This scalar product is linear in the variable ip and antilinear in the variable ip. 
Throughout this exposition we assume that d is finite. Denote by B(H) the ★-algebra of 
all operators on H. The complex <i 2 -dimensional vector space B(H) will also be viewed 
as a Hilbert space with the scalar product (X\Y) = TrX+y where denotes the 
adjoint of the operator X. Denote by S(H) C B(H) the compact convex set of all 
nonnegative (definite) operators of unit trace. Any element p in S(TC) is called a state 
of the system. The extreme points of S(7i) are precisely one dimensional orthogonal 
projections. They are called pure states. In the Dirac notation any pure state can be 
expressed as \ip >< ip\ where if) is a unit vector in 7i. Denote by V{7i) the set of all 
orthogonal projection operators (or, simply, projections) on 7i. Any element P in V(7i) 
is called an event concerning the system and the quantity TrpP is interpreted as the 
probability of the event P in the state p. In the context of quantum information theory 
the state of a quantum system can be utilized as an information resource. If the system 
is in an unknown state p it is important to estimate p from "independent repeated 
measurements". If we choose and fix an orthonormal basis {eo, e±, . . . , ej-i} in 7i then 
p is described in this basis by a nonnegative definite matrix ((pij)) where p^ = {ei\p\ej). 

Thus determination of p involves the determination of d 2 — 1 real parameters, namely, 

d-i 

pa, i — 1, 2, . . . , d — 1, Repij, Imp^, 0<i<j<d— 1. (Note that p o = 1 — Yl Pa an d 

i=i 

Pij = Pji-) 

By an elementary measurement M. = {Po, Pi, ... , Pd-i} we mean a family of d mu- 
tually orthogonal one dimensional projection operators Pj, j — 0, 1, 2, . . . , d — 1 so that 

d-i 

Y Pj = I, the identity operator. If the measurement Ai is performed when the state 
o 

of the system is p, the result of such a measurement is one of the classical outcomes 
j G {0, 1, 2, . . . , d — 1} with probability TrpPj = pj for each j. Independent repeated 
trials of the measurement in the same state p yield frequencies fj for each elementary 
outcome j and fj can be viewed as an estimate of pj for each j. Thus an elementary 

measurement covers at most d— 1 degrees of freedom concerning p in view of the relation 

d-i 

^2 pj — 1. In order to estimate p it is therefore necessary to examine the frequencies of 

j=o 

elementary outcomes in at least d + 1 elementary measurements A4j, < j < d where 
no two of the measurements M.i and M.j have any "overlap of information". Such an 
attempt is likely to cover all (d+ l)(d— 1) = d 2 — 1 degrees of freedom involved in recon- 
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structing or estimating the unknown p. To bring clarity to the notion of "nonoverlap of 
information" in a pair of elementary measurements it is useful to look at the ★-abelian 
algebra 



A(M) 




cij e C,j — 0, 1, . . . , d — 1 



d-i 

Any element X — J] XjPj in A{M) can be looked upon as a complex- valued observable 
j=o 

where Pj is interpreted as the event that "X assumes the value . Of course, this 
is justified if all the distinct scalars. If x is any scalar then the event that X 

assumes the value x is the projection Yl Pj- Thus the subalgebra CI C A(M) consists 

J'.Xj=X 

precisely of constant-valued observables. Such an interpretation motivates the following 
formal definition. 

Definition 1.1 Two elementary measurements M = {P , Pi, ■ ■ ■ , Pd-i}, 

■M.' = {Qo, Qi, ■ ■ ■ , Qd-i} are said to be weakly mutually unbiased (WMUB) if 

A(M) n A(M') = CI, 

and strongly mutually unbiased (SMUB) if, in the Hilbert space B(H), the subspaces 
A(M) CI and A(M') CI are mutually orthogonal. (Here, for two subspaces £1 C 
S2 C B(H), S2 S\ denotes the orthogonal complement of Si in S2) ■ 

Clearly SMUB implies WMUB. We shall now describe these two properties in terms 
of the quantities TrPjQj. 

Proposition 1.2 Two elementary measurements M = {P , Pi, ... , Pd-i}, M.' = {Qo, Qi, ■ ■ ■ , Qd-i} 
are SMUB if and only if 

Tr PiQj = (T 1 for all i,j,E {0, 1, 2, . . . , d - 1}. (1.1) 



Proof: Note that the subspaces A(M) CI and A(M') CI are respectively spanned 
by the subsets {Pj - d' l I, < j < d - 1} and {Qj - d'H, < j < d - 1}. Thus the 
orthogonality of these two subspaces is equivalent to the condition 

0=(P i - d- l l\ Qj - d- l I) = Tr(P; - - (T 1 /) = (TrPiQj) - (T 1 
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for all i,j in {0,1,2,..., d- 1}. □ 

Proposition 1.3 Let M = {P , Pi, ... , Pd-i}, M! = {Qo, Q±, . . . , Qd-i} be two elemen- 
tary measurements. Suppose 

L = [Tr(P t - P )(Qj - Qo)} , i, j e {1, 2, . . . , d - 1} 

and J d -i is the (d — 1) x (d — 1) matrix all the entries of which are unity. Then M and 
M are WMUB if and only if 

det + J d _i + d^LJa^l) - W) > 0. (1.2) 

Proof: Let X e A(M) fl „4(.A/f )• Then there exist scalars a,, fy, i, j e {1, 2, . . . , e? - 1} 
such that 

d-i 

X = d-\TrX)I + Y, a ^ P i- p ») 

i=\ 

d-l 

= d-\TrX)I + h AQj ~ Qo)- 

3=1 

Thus M and M' are WMUB if and only if the set {P 1 - P , P 2 - P , . . . , P d - P , Qi - 
Qo,Q2 — Qo,---,Qd — Qo} of 2(d — 1) elements in the Hilbert space B(7i) is linearly inde- 
pendent. This, in turn, is equivalent to the strict positive definiteness of the partitioned 
matrix 



" [Tr(P, 


-Po)(Pj 


~Po)] 


[Tr(P, 


-Po)(Qj 


- Qo)] ' 


. [Tr(Qi 


-Qo)(Pj 


-p )] 


[Tr(Q, 


-QoKQi 


- Qo)] . 



,i,je{l,2,...,d-l} 



of order 2(d — 1). We have 

Tr(P - P Q ){Pj - P ) = Tr (Qi - Q ){Q j - Q ) 
Thus, M and M' are WMUB if and only if 



Id-l + Jd-1 


L 




Id-i + Jd-i 
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has a strictly positive determinant. Left multiplication of this matrix by the matrix 



' Id-1 


-L(I d _i + Jd-i)- 1 ' 





Id-i 



with unit determinant yields the equivalent condition 

det + J d _! - L(/ d _! + J^y 1 ^) > 0. (1.3) 

Since 

(Id-1 + ^d-l) 1 — Id-1 ~ d Jd-li 

condition ()1.3|) reduces to condition ()1.2|) . □ 

Corollary 1.4 If the matrix L of Proposition 1.3 satisfies the inequality ||L|| < l(where 
|.| is the standard operator norm in the ^-algebra ^(C^ 1 ) then M, M' are WMUB. 
Furthermore M., M.' are SMUB if and only if L = 0. 

Proof: Immediate. □ 

In the context of minimizing the number of elementary measurements required for 
estimating the state p of a quantum system Proposition 1.2 emphasizes the importance 
of the search for d + 1 elementary measurements which are pairwise SMUB. When d is 
a prime power p r the existence of such a family of SMUB measurements was proved by 
Wootters and Fields [7]. Alternative proofs of this result were given by S. Bandyopadhyay 
et al in [3] and Pittenger and Rubin in [6]. In this paper we shall present a proof of 
the same result by using the commutation relations of Weyl operators in the L 2 space 
of the finite field Fpr. When d = p^p™ 2 ■ ■ - p™ n with p^s being prime we shall use the 
Weyl commutation relations in the L 2 space of the additive abelian group ®f = i^ Pi ^ and 
study the problem of estimating the unknown state of a <i-level system. This leads to an 
interesting reconstruction formula for a state p in terms of probabilities of d 2 — 1 events 
arising from YTi=i (pT* + -0 elementary measurements. However, one would like to express 
p in terms of the probabilities of elementary outcomes in (nr=i + -0 elementary 
measurements. 
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2 The case d = p r 



Let dim Tt = d = p r be a prime power. For any prime power q denote by ¥ q the unique 
(upto a field isomorphism) finite field of cardinality q. Choose and fix any nontrivial 
character \ of the additive group ¥ d and put 

(x,y) = xN> x,y G ¥ d . (2.1) 

One can, for example, look upon ¥ d as an r-dimensional vector space over ¥ p , express 
any element x in F^ as an ordered r-tuple: x — (si, s 2 , . . . , s r ) where < Sj < p — 1 for 
each i and put 

X(V)=exp si. (2.2) 

P 

Then we have | <cc, y> | = 1, = {x,y x + y 2 ) = (x,yi)(x,y 2 ) and a; = if 

(x, y) = 1 for all y in F^. In other words, (., .) is a nondegenerate symmetric bicharacter 
for ¥ d . Identify the Hilbert space Ti with L 2 (F^), using the counting measure in F^, and 
put 

\x >= l {x} ,xe ¥ d 

where 1{ X } is the indicator function of the singleton subset {x} in ¥ d . Then {\x >, x G ¥ d } 
is an ortho normal basis for H labelled by the elements of ¥ d . Now, consider the unique 
unitary operators U a , Ub in 7~C determined by the relations 

U a \x > = \a + x >, 

V b \x > = < b, x > \x > for all x £ ¥ d . 

Then we have 

U a U b = U a+b ,V a V b = V a+b , (2.3) 

V b U a = (a,b)U a V b . (2.4) 

Elementary algebra shows that 

(2.5) 

for all di, a,2, bi, b 2 in ¥ d . In particular, the family {U a V b , a, b G F^} of d 2 unitary operators 
constitute an orthogonal basis for the Hilbert space B(H). This is an example of a unitary 
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error basis in the theory of error correcting quantum codes [4] . Notice also the fact that 
{U a } and {H} are like the position and momentum representations obeying the Weyl 
commutation relations in classical quantum mechanics. In view of this property we call 
any operator of the form XU a Vb, \X\ = 1, a, b G F d a Weyl operator. We say that (|2.3|) 
and (j2.4j) constitute the Weyl commutation relations. The usefulness of such an error 
basis of Weyl operators in the study of quantum codes has been explored in [1], [2], [5]. 
We shall slightly modify the error basis {U a Vb} by multiplying each element U a Vb by an 
appropriate phase factor. Once again viewing as an r-dimensional vector space over 
F p , expressing any x G F d as an ordered r-tuple x = (si, s 2 , . . . , s r ) with < Sj < p — 1 
for each i and considering the basis elements = (0, 0, . . . , 0, 1, 0, . . . , 0) of the field F^ 
with 1 in the i-th position and elsewhere we write x = s^i + ^2^2 + • • • + s r e r and 
define 

a(a,x) = x ^" -i + 5Z~^~2 — ' ^ ^ ¥d 

where x is the character chosen and fixed at the beginning of this section. 
Now put Frf = Fd U {00} and write 

„ TI v ( a{a,x)U x V ax if a G F d , x G F d , 
W{a,x) = < (2.7) 
I V x if a = 00. 

Then we have the following proposition. 

Proposition 2.1 The family |/, W(a, x),a G F^, x G F^ \ {0}} is an orthogonal basis of 
unitary operators for the operator Hilbert space B(TC) satisfying the relations 

W{a, x)W{a, y) = W{a, x + y) for all ae¥ d ,xe F d . (2.8) 



Proof : The first part is immediate from the fact that the family of operators under 
consideration differs from the family {U x V y , x, y G F^} only by a scalar factor of modulus 
unity in each element. If a G F d , x = s i e i, V — Yl ^i e i we have from (|2.3|) (|2.4|) 

W(a, x)W(a, y) 
= a(a, x)a(a, y) (ax, y)U x+y V a ^ x+y) 
= a(a, x)ac(a, y)a(a, x + y)(ax, y)W(a, x + y) 
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where the coefficient of W(a, x + y) is of the form x{ az ) with 

z = SiSjeiej + — - ej + titjeiej + — - 

- J2( Si + U)(s, + t,)e,e, - £ fo + ^+'j- 1 ^ + £ *f 
= 0. 

This proves when a G F d . When a = oo, (EH) is a part of (|2~3)l . □ 

Theorem 2.2 There exists a family of one dimensional orthogonal projection operators 
{P(a, x),a G Fd, x G F^} satisfying the following : 

(i) W(a,x)= £ (x,y)P(a,y) 



(ii) P(a,y) = d- 1 £ (z.y^a,*), 

(iii) P(a,x)P(a,y) = 6 x , y P(a,x), 

( iv ) ExeF d &) = J > 

(v) Tr P(a, x)P{b, y) = cT 1 for all a ^ b; a, b G F d ; x, y G F d . 

Proof : By Proposition 2.1 the correspondence x — > W(a, x) is a unitary representation 
of the additive abelian group F d and {(.,y),y G F d } is the set of all its characters. 
Thus the decomposition of {VF(a, .)} into its irreducible components yields a spectral 
measure P(a, .) on F^ satisfying (i), (iii) and (iv). Substituting from (i) the expression for 
W(a, x) in the right hand side of (ii) and using the orthogonality relations for characters 
we get (ii). Taking trace on both the sides of (ii) and observing that W(a, 0) = I and 
Tr W(a,x) = for x ^ we get Tr P(a,y) = 1. Thus each P(a,y) is a one dimensional 
projection. Substituting for P(a, x) and P(a, y) from (ii) in the left hand side of (v) we 
have from lf277|) . (1231) and (12~3J) 

Tr P(a,x)P{b,y) 
= d~ 2 (x,z 1 )(y,z 2 )TrW(a,z 1 )W(b,z 2 ) 

Zl,Z 2 & d 

= cT 2 ( x ' z i)(yi z 2) a ( a i z i) a & z 2 )(az 1 ,z 2 )TrU zl+Z2 V azi+bz2 . 

Zl,Z2&d 
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Now observe that the (zi, z 2 )-th term of the sum on the right hand side is nonzero only 
if z\ + z 2 = 0, az 1 + bz 2 = 0. If a ^ b this is possible only if z± — z 2 — 0. This proves (v). 
□ 

Corollary 2.3 Let M a = {P(a,x),x G ¥ d }. Then {M a ,a G ¥ d } is a set of (d + 1) 
elementary measurements which are pairwise SMUB. 

Proof: Immediate from Propostion 1.2. □ 

Our next result yields a recovery formula for any state p from the probability distri- 
butions {Tr pP(a, x),x G F^} on ¥ d arising from the measurements {A4 a , a G F^}. 

Theorem 2.4 Let {P(a : x),a e¥ d ,xE ¥ d } be the projections in Theorem 2.2. Then, 
for any state p on L 2 (¥ d ) the following holds: 

(i) P= E Ef&^Ml-^M) 

(ii) p= E (^){TrpF(a, ?/ )}Vy(a, a :) 

Proof: From the first part of Proposition 2.1, it follows that p admits the expansion 
p = d- 1 \l+ [TrpW(a,xy] W(a,x) 

aew d 
K xev d \{o} 

in terms of the orthogonal basis arising from the Weyl operators. Now substitute in the 
right hand side the expressions for W(a, x) in (i) of Theorem 2.2 and use the orthogonality 
relations for characters: 

Then we obtain the identity (i) of the theorem. If we substitute for P(a, z) from the 
identity (ii) of Theorem 2.2 we obtain the second identity of the theorem. □ 

Remark: If we make repeated independent measurements A4 a , obtain the frequencies for 
the different events P(a, z) and substitute those frequencies for the different probabilities 
Tr pP(a, z) in the unknown state p we will get an unbiased and asymptotically consistent 
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estimate p of p but p may not be a positive operator. One may replace p by the normalised 
version of the positive part or the modulus of p at the cost of losing unbiasedness. This 
also increases the computational cost. 

3 Estimation of states in the general case 

Let d = p^p™ 2 ■ ■ ■ Pn n be the decomposition of d into its prime factors p% < p<i < 
■■■ < p n . Write dj = p™ 3 . We may identify the d- dimensional Hilbert space Ti with 
Ti\ <E> Ti 2 <8> • • • ® Ti n where Tij = L 2 (F dj ), ¥ d . being the finite field of cardinality dj. 
Following the definition in (J2.7|) construct the unitary operators W^\aj, Xj) when d = 
dj, j = 1,2, ...,n and using Theorem 2.2, the corresponding projections P^\aj,Xj), 
where dj G F<j., Xj G F^.. We now adopt the following convention: for any operator 
X in L 2 (¥ dj ) = Hj denote by the same symbol X the operator in Ti. defined by X = 
X\ ® X2 <E> • • • <E> X n where Xi is the identity operator in Tii when i 7^ j and Xj = X. 
The operator X thus defined in Tt = Hi ® <E> • • • <8> Tin is called the ampliation of X 
in ft,- to H. Since can be identified with B(Hi) ® fi(W 2 ) ® • • • ® B(W„) as Hilbert 

spaces as well as *- algebras it follows from Proposition 2.1 that the family 

T = {l,W^\a il ,x il )W^\a ia ,x ia )---W^\a ir ,x ir ), 
<Hi e F dj . , x< . G F d . \ {0}, j = 1, 2, . . . , r, 

1 < ii < i 2 < • • • < i r < n, r — 1, 2, . . . , n} (3.1) 

of unitary operators in Ti constitute an orthogonal basis for the operator Hilbert spaces 
B(TL). Note that the cardinality of T is, indeed, equal to 

n 

1+ E E 

r=l l<ii<i2<-"<ir<n 

= (i + d?-i)(i + 4-i)... (1 + 4-1) 
= «...d 2 n 

= rf 2 , 



10 



the dimension of B(TL). For any subset J = {11,12, ■ ■ ■ , v} C {1, 2, . . . , n} where 1 < i\ < 
%i < ■ ■ ■ < i r < n, define 

(i(t/) d^d^ . . . &{ r 

d'(j) = k + + + 

and for any state p in TC, put 

s p( j ) = E {^pP {n \^y l ,)P {l2 \^y l2 )---P {ir \a lr ,y ir )} 

3 l 3 
3 % 3 

P^\a^y n )P^\a l2 ,y l2 ) . . . P^(a ir , y ir ) (3.2) 

where {P^(aj, yi)} are the one dimensional projections in Tii determined by the unitary 
representation Xi — > W^'(ai, Xi) of the additive group F^ according to Theorem 2.2 and 
ampliated to the product Hilbert space 7i = l~ii®TL2®- ■ •®'H n - Thus S P (J) is an operator 
in H determined by the probabilities Tr pP^ 1 ' (a^y^) P^ 2 \ai 2 yi 2 ) . . . P^ r \a ir y ir ) and the 
projections P^^a^, t/ijP^^a^, y i2 ) . . . P M (a ir ,y ir ) of dimension rL,^^...^}^ with 
aj's varying in F^ and y^s in F^ for any i. With these notations and the convention 
S p (4>) = I, we have the following theorem for the recovery of p from the probabilities. 

Theorem 3.1 Let p be any state in 7i. Then 

P = Yl (-ir |J| *VJ) (3.3) 

JC{l,2,...,n} 

where S P (J) is given by ()3.2|) and |J| is the cardinality of J. 

Proof: Since the family JF of unitary operators in (J3.1j) is an orthogonal basis for B{7i) 
we can expand the state p in this basis as 

p = {d 1 d t ...d n )- 1 \l + Y i £ £ 

r=l l<i!<i2<-<V<n 0i e v d ,Xi.eW d \{0} 

[Tr pW^\a h ,x n y ■ ■ - W^^^x^} W^\a h , x h ) ■ ■ ■ W^\a ir , x ir )} .(3.4) 
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From Theorem 2.2 we have for any fixed i 



J2 W®(a i ,x i y®W®(a i ,x i ) 

Xi& di \{0} 

Jx~y)(x i ,z)P^(a i ,y)^P^(a i ,z) 

»=i6F d .\{0} 

= d t P^Kv) ® P {l) (a iiy ) - /« (g) /«, 

being the identity operator in Tti. Using this identity and elementary properties of 
relative trace, equation (|3.4j) can be written as 



P = (d 1 d 2 ...d n )- 1 ^^(-l)W-Wd{K)d!{J\K) 

J KcJ 

x J2 i Tr P pikl } K > ^ ) ^ ( a *> > ^ ) • • • » W.) } 



xP^\a kl ,y kl )P^\a k2 ,y k2 ) ■ ■ ■ P {K \a ks , 



where J varies over all subsets %\ < i 2 < ■ ■ ■ < i r of {1, 2, ... ,n} and K varies over all 
subsets k\ < k 2 < ■ • ■ < k s of J. Now using the definition in (|3.2|) we can express p as 

p = < K )S P {K) 

KC{l,2,...,n} 

where 

a{K) = (d 1 d 2 ...d n )- 1 d(K) (- 1 )' L|rf '( L ) 

L:LnK=<t> 

= (_i)«-l A 'l.n 

Remark From Theorem 3.1 it is clear that p is recovered from the probabilities for the 
elementary events 

p( 1 \a 1 ,x l )p( 2 \a 2 ,x 2 )...P {n) (a n ,x n ), a, e W dt Xi e ¥ di . 

In other words the determination of p involves (d% + l)(d 2 + 1) • • • (d n + 1) elementary 
measurements. As mentioned in the introduction one would like to determine p by 
d\d 2 . . . d n + 1 measurements. 
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